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Abstract. Moving vortex matter, driven by transport currents independent of time, in which 
vortices and anti- vortices coexist is investigated theoretically in thin superconducting films with 
nanostructured defects. A simple London model is proposed for the vortex dynamics in films 
with periodic arrays of nanomagnets or cylindrical holes (antidots). Common to these films is 
that vortex anti- vortex pairs may be created in the vicinity of the defects by relatively small 
transport currents, because it adds to the current generated by the defects - the nanomagnets 
screening current, or the antidots backflow current - and may exceed locally the critical value for 
vortex anti- vortex pair creation. The model assumes that vortex matter dynamics is governed 
by Langevin equations, modified to account for creation and annihilation of vortex anti- vortex 
pairs. For pair creation, it is assumed that whenever the total current at some location exceeds 
a critical value, equal to that needed to separate a vortex from an anti-vortex by a vortex 
core diameter, a pair is created instantaneously around this location. Pair annihilation occurs 
by vortex anti- vortex collisions. The model is applied to films at zero external magnetic field 
and low temperatures. It is found that several moving vortex matter steady-states with equal 
numbers of vortices and anti- vortices are possible. 



The dynamics of vortex matter in thin superconducting films with nanostructured defects has 
attracted a great deal of attention recently [T|. Most works done so far assume that the vortex 
matter is made up of identical vortices with unit vorticity, and that the laws governing its 
motion are those for ordinary matter consisting of identical classical particles with zero mass 
[2j. However, an important difference between vortex matter and ordinary matter is that 
vortex anti-vortex (v-av) pairs can be created by thermal fluctuations or by supercurrents, and 
annihilated by collisions [3J. These processes can be neglected only at low temperatures, and for 
supercurrents everywhere small compared to the critical current for v-av pair creation. These 
conditions are well satisfied in homogeneous films at low temperatures, and when the transport 
current is uniform, since it cannot exceed the critical current, otherwise superconductivity is 
destroyed by proliferation of v-av pairs. However, in films with nanostructured defects, the 
super current is not in general uniform, and may exceed the critical current locally, creating only 
a few v-av pairs. As a consequence, superconductivity is not destroyed, and states of moving 
vortex matter with coexisting vortices and anti-vortices are produced. These states exist even in 
the absence of applied magnetic fields. Two systems are considered here: i) films with arrays of 
nanomagnets, for which the transport current is uniform, but a non-uniform screening current 
flows, induced by the nanomagnets. ii) films with arrays of holes, for which the transport 
current itself is non-uniform, because it cannot go through the holes. This paper demonstrates 
the existence of such states at zero applied fields, and studies some of its properties. 



The theoretical model is the following. Vortex matter in thin superconducting films, 
characterized by thickness d, penetration depth X ^ and vortex core radius ^, interact 
with the following defect arrays: i) a square lattice of magnetic dipoles, with magnetic moment 
m parallel to one of the unit cell sides (m = mx ) , placed at a distance zq from the film surface 
( FigH]). ii) a square lattice of cyhndrical holes (antidots) of radius R ^ { FigH]). In both 
cases the defect lattice unit cell has dimensions x a^. The vortex matter is modeled as a 
system of classical particles with zero mass. To account for v-av pair creation and annihilation, 
it is assumed that the time evolution can be broken into two parts: motion of the vortices and 
anti-vortices, and creation of v-av pairs. The first part is governed by Langevin equations ^ 

V^-qj{F] + F'^ + Fj)+Tj, (1) 

where 77 is the friction coefficient, j labels vortices and anti- vortices present at time t, Tj{t) and 
Qj = ibl,ili2, ... denote, respectively, the corresponding position vectors and vorticities (q > 
for vortices, g < for anti- vortices) . FJ, and Fj are, respectively, the forces of interaction of 
a test vortex (q = 1) located at rj{t) with the other vortices and anti- vortices, and with the 
defect array. Fj is the Lorentz force of the transport current, J, on a test vortex at rj(t), that 
is Fj = {(f)od/c)3 X z. Tj is the random force. In the film with the dipole array, J is uniform 
and makes an angle a with the ?/-axis ( FigH]). In both cases J is independent of time. The 
forces and F^ are also the Lorentz forces on a test vortex due, respectively, to the currents 
generated by the other vortices and anti- vortices, J^(r), and by the defect array, J(^(r). These 
equations naturally account for v-av pair annihilation, which happen when a vortex and an 
anti- vortex collide. For v-av pair creation, the model assumes that if, at a given instant of time, 
the current at some location exceeds a critical value, a v-av pair is added around this location 
instantaneously. 

The force F^ is given by FJ = —^j^iqiUvv{rj — r^), where Uvv{t) is the vortex- vortex 
interaction potential in the London limit. Essentially, Uyy{r) ~ — eodln(r/A) (eo = ((/)o/47rA)^, 
and A = 2A^/(i), for r > ^. For r < ^, a Gaussian cutoff is used in Uyy^ in order to account 
for the vortex core [4J. Creation of a v-av pair by a supercurrent occurs if its Lorentz force 
can separate a vortex from an anti- vortex by a distance ^ 2^, along the direction of the 
Lorentz force. The critical current, J^, is estimated by balancing the energy change for pair 
creation, AE = 2£^sj + 2Uyy{2^) {Egf ^ eodln(^/A) = vortex self energy), with the work 
done by the Lorentz force, 2^(j)od Jj^/c. The result is ^ (ceo/(/>oO ln2 = 0.9 Jd, where 
Jd = C(/)o/(12v^7r^A^^) is the depairing current. The relevant current for v-av pair creation is 
the total current, Jt = J + + J^;- When 3t exceeds locally, v-av pairs are created. The 
detailed procedure for v-av pair creation is discussed next for each defect array. 
Dipole Array. The vortex-defect force is the magnetic interaction between the vortex and the 
dipoles, given by F^(r) = — VX]j^^ZY(r — R^) where U{r) = (egd) (47rm/(/)o^o) (^o^^ • x/r^)(l — 
(1 + r^/zg)"^^^) is the potential of a single dipole, r is the vortex position vector with respect 
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Figure 1. Film and defect arrays: a) Side view of dipole array, b) Top view of dipole array, 
and transport current orientation, c) Top view of antidot array. 



to the dipole position (projected in the film plane), and denotes the dipole array translation 
vectors. The screening current follows from = {(j)od/c)3d x z. The resulting for a dipole 
array with » zq is shown in Fig|2^, close to the position of one of the dipoles, where is the 
largest. The following procedure is used to add v-av pairs. If at a point rg, Jr(ro) exceeds 
(Jt(^o) > Jm)^ a v-av pair is added around rg: a vortex = 1) at rg + 5, and an anti- vortex 
at ro — 5, where 5 has modulus ^ and is oriented parallel to the Lorentz force of JT(ro), that is 
5 — ^JT(ro) X i/ Jt{yo). A fundamental assumption of the model is that v-av pairs can he added 
in the manner described above only around points fq for which \ J + J(i(ro) |> Jm- These points 
are localized in regions close to the dipole positions, or between them, depending the orientation 
and magnitude of J. An example is shown in Figl2}3. This assumption is necessary in order to 
avoid unphysical behavior caused by the vortex core. What happens is that the supercurrent 
generated by a vortex or an anti-vortex has values close to in the proximity of the core. 
When this current adds to another - transport, screening, or from another vortex or anti- vortex 
- the resultant may exceed in the proximity of the core, leading to creation of v-av pairs. 
This is unphysical. What should happen in a more realistic model is deformation of the vortex 
core. A more detailed justification of this assumption will be given elsewhere. 
Antidot Array. Here only a single antidot is considered in detail, which is relevant for arrays 
where the antidots are far apart (a^ >> R ^ A vortex, with vorticity interacts with the 
antidot through the two images it induces on the hole: one, with vorticity at the center of the 
hole (r = 0), the other, with vorticity — at r/ = (i?/r)^ r, where r is the vortex position vector 
0. Here this interaction is included in F^, since it depends on the vortex matter distribution. 
Thus, 3y also contains the current generated by the vortex images outside the holes, and F^ = 0. 
The transport current is the sum of an uniform current Jq, with the backflow current J^j, that is 
J = Jo + J5/. For the circular hole, J^j has a dipolar spatial dependence, identical to that of an 
ideal fluid. An example of J is shown in Fig|2]3. The maximum values of J occurs at two points 
located in the circumference of the hole perpendicular to Jq. At these points J = 2Jo. Thus, 
when Jo >= O.SJ^, J exceeds in regions located in the vicinity of these points, and v-av 
pairs are created [7J. A vortex (anti- vortex) is added in the part of the region where J > 
and the Lorentz force on the vortex (anti-vortex) points away from the hole center. An example 
of these regions is shown in Figl2]i. The critical current for for v-av pair creation is, essentially, 
Jjn in the limit ^. Thus, the final result of this process is creation of a v-av pair near the 
defect, as in the film with the dipole array. 

Numerical solutions of Eq. ([T]) are obtained at low temperatures for the film with the dipole 
array, using a discrete version of Eq. ([T]) on a square grid. The results confirm the existence of 
a rich variety of steady-states with equal numbers of vortices and anti- vortices. Typical ones 
are shown in FigjSl for J < J^, 90^ < a < 180^, and m chosen so that the screening current 
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Figure 2. a) Screening current around a dipole (arrow), b) Regions below a dipole where v-av 
pairs are created (red), for a = 135^, J/Jm — 0.55. c) Transport current in the vicinity of a 
circular hole of radius R for Jq in the negative ^/-direction, d) Regions near the hole where a 
vortex (red +) or anti- vortex ( red — ) are created, for J/Jm — 0.75. 




Figure 3. Results for a film with dipole array (m = O.SS^oC, = 2^, = 16^)- a) Time 
averaged number of vortices per dipole vs. transport current, b) Number of vortices per dipole 
vs. time (arbitrary units), c) Trajectories of vortices (light) and anti- vortices (dark). Circles: 
regions where v-av pairs are created (VC) and annihilated (VA). 

is everywhere smaller that J^. The number of v-av pairs per dipole, Np, averaged over time 
(denoted < Np >) depends on J and a as shown in FigjStt. The time dependence of Np, shown 
in FigjSb, is oscillatory, varying between one and two, almost periodically in this particular 
case. Typical trajectories for the vortices and anti- vortices and the regions where v-av pairs 
are created and annihilated are shown in FigjSt. For a < 80^ no steady-states are found for 
J < Jjn [8j. The behavior shown in Figl3]is a direct consequence of the model basic assumptions. 
The dependence of < Np > on a results from the spatial dependence of the screening current 
(Fig|2^): for a > 90^ the transport current and the screening current add near dipoles, favoring 
v-av pair creation, whereas for a < 90^ the opposite occurs. The time-dependence of Np is 
due to oscillations in time of J^, caused by the motion of the vortex matter. During the time 
interval when Np(t) = 1, < Jm in the regions where v-av pairs are created, but increases 
in these regions as the vortices and the anti- vortices move away from the dipoles. Eventually 
Jt > Jm in parts of these regions and one new v-av pairs is created by each dipole. After that 
Np(t) = 2, until collisions between vortices and anti- vortices, taking place between the dipoles 
(FigjSt), reduce it to one. 

Steady-states of moving vortex matter with similar properties are expected to exist in films 
with periodic arrays of holes, as well as in films with non-periodic defect arrays. 
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